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Monte Carlo data simulating phase transitions in Ising strips D x L, (D <C L) with periodic 
boundary conditions show that T C (D) = for D < D* ~ 6 and < T C (D) < T c (d = 2) for D > D*. 
Regular scaling of ML 13 '' 1 ' vs \T — T C \L}^ V is obtained only for D > D*and the Monte Carlo effective 
susceptibility critical exponent 7e//(-D) is shown to be well described by 7(d) = /3(d) [5(d) — 1] with 
d eff (D) given by d eff (D) ~ 1.5 + (^)(D - 6) and (3(d) = (ff - f), ^(d) = (ff - I), which can 
be understood as valid with d e ff(D). 



PACS numbers: 64.60.-i, 68.18.Jk, 64.60.Cn. 

I. INTRODUCTION 

It is well known [l| that systems at criticality, for in- 
stance Ising magnets at (T — ► T C1 H — > 0), are usually 
not only scale invariant but also conformally invariant, 
i.e. the pertinent local transformations look like combi- 
nations of dilations, rotations and traslations, but shear 
distortions are not allowed. This fact has been used to 
prove that in two dimensions, the Ising critical exponents 
must be fractional. 

Monte Carlo methods Q have been extensively used 
to characterize phase transitions in Ising systems and are 
used in this work to stidy D x L strips. 

Recent Monte Carlo |3j data on large three dimen- 
sional Ising lattices (L 3 ~ 1.5 x 10 6 spins) appear to 
indicate that for 2 < d < 4 (including d — 3) the critical 
exponents are either well defined fractional Q or indis- 
tinguishable from fractional, and that they are determi- 
nated by compact expressions for j3(d) and 8~ 1 (d). (See 
Figur(0J 

In the present work we carry out an empirical investi- 
gation of the critical behaviour of strips D x L (with 
increasing width D <C L) by means of Monte Carlo 
simulations of the phase transition for long strips with 
D = 2,3,4... up to D = 52, using periodic bound- 
ary conditions, to ascertain whether the effective critical 
exponents (3(d) and 5~ x (d) describing the phase transi- 
tion are given by the same expressions linear in d which 
describe well the behaviour for higher dimensionalities 
(d = 2,3,4). 



II. MONTE CARLO RESULTS 

Numerical finite size simulations of the phase tran- 
sitions in D X L Ising systems with periodic bound- 
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FIG. 1: Critical exponents (3(d) = (ff - \) and 5~ 1 (d) = 
(ff — |) extrapolated from the respective fractional values 
[(3(4) = 1/2, (3(2) = 1/8] and [5~\4) = 1/3, S~\2) = 1/15]. 



ary conditions were performed for D = 2 to D = 52. 
Wolff cluster algorithms were used in strips of length 
250 < L < 5000. Periodic boundary conditions were 
used always in the L direction but they had not much 
effect for large L values as it is to be expected. The ther- 
malization time, relaxation time and number of states 
were increased steadily until the final results were not 
appreciably affected by further increases. 140000 Monte 
Carlo steps per spin for each temperature were taken. To 
reduce the critical slowing down at T ~ T c a single clus- 
ter Wolff algorithm was used. Initial conditions at a any 
given temperature in a closely spaced set of temperatures 
(AT < 0.01) were taken from the equilibrium conditions 
at the previous temperature. 

The critical temperature of strips with width D <C 
L which is a non-universal quantity, was determined by 
means of the Binder Cumulant method 0] , crossing data 
for D x L with data for D x 2L. It was found that for 
D < 6 no crossing took place at T > 0, and it was checked 
whether scaling of ML^I V vs ^[L 1 ^, with |e| = \T-T C \, 



2 



was possible using (3 = and v < 1/2. It was found that 
no such scaling took place, but, what could be called 
"one — dimensional" scaling 1], with |e| = e _4 / T was 
observed to take place for D = 2,3, 4, 5. 

For D > 6 the Binder Cumulant method did provide 
non-cero transition temperatures given by [2| 



D* ~ 6, Equation Q gives 7(d) in terms of /3(d) and 
S~ 1 (d), as specified by Equations J3J) and As shown, 
j e ff(D) is given by 



l*ff{D) = G 



A 



D-D" 



— B (D - D*) , D > D* (6) 



T C (D) = T c (d = 2)[l-e 



(1) 



with D* ~ 6, and to ~ 0.353 ± 0.011; and regular 
scaling of Mtf'" vs leli 1 ^, with |e| = \T - T c \ was 
observed to hold, with (3(d) and v(d) evolving smoothly 
between the respective values for d = 1.5 and d ~ 2 as 
specified below. 

To describe the evolution of the critical exponents as a 
function of two-dimensional strip width D we note first 
that f3{d), S-\d) and 7(d) = (3(d) [6(d) - 1] for d = 2, 
and d = 4 are given by fractional values specified by 



5- x (d) 



2d 
15 



7 (d) = ^(d)[«y(d)-i] = y+^ 



2d -3 



3 

16 



(2) 



(3) 



2d - 3 



(4) 

For d = 3 Monte Carlo data Q are consistent within 
narrow statistical error bars with fractional exponents 
(3(d = 3) = ^, ^(d = 3) = § and 7 (d = 3) = §, as 
given by Equations J2J> and (@J respectively. 

Figure shows plots of x" 1 vs |e| = |e~ 4 / T | for 
D < D* = 6 and Figure [3] shows like plots of x 1 vs 
|e| = \T - T c \ for D > D*. It can be seen that j e ff(D) 
diverges at a width £> ~ 6 from both sides. From these 
data the evolution of 7e / / (D) can be determinated with 
fair accuracy, and it can be used in conjunction with 
Equation J3J to determine the effective dimensionality of 
strips with growing D. 

This is done in Figure 21 where a plot of d e / / (D) re- 
sults in a quiasi-linear dependence of d e ff(D) starting 
at d e ff(D*) — 1.5 (corresponding to 6~ 1 (d) = af- 
ter Equation |2J and growing up from this value towards 
d e f f (L) = 2, i.e. for rectangular LxD strips with D — > L 
(approaching square lattices). d e ff(D) therefore, is ap- 
proximately given by 



(5) 



d eff (D) = 1.5+ [ — ) (D-D* 



where D* ~ 6. 

Figure gives a plot of j e ff(D) which shows directly 
how this critical exponent blows up as D approaches 



with G ~ 11/8, A ~ 1125/24 and B ~ 3/3200, in 
good agreement with 7(d) given by Equation For 
D < D*, i.e. for D = 2,3,4,5, 7e//(-D) grows up steeply 
with D and blows up at D ~ 6, which corresponds to 
(3(d c ) = where d c = 1.333..., as given by Equation 
7e//(D) for D < D* , i.e. for the interval 1 < d < d c 
where 7 ~ (3 ■ 6 ~ (0) • (00) , is undefined, but it can be 
fitted reasonably well by 



j e ff(D) ~C-(D-D*)-\ D < D* 



(7) 



with C ~ 4, which is compatible with 7 (d = 1) = 1/2. 

The hyperscaling relationships [| result in a critical 
exponent v describing the temperature dependence of the 
correlation length £ which is related to the dimensionality 
d and to other exponents by 



V -\d)= Vl (d) = 



f3(d)[5(d)-l] 



(8) 



For D > D* the denominator in the right hand side 
of Equation JSJ can be written as {7(d) + 2 /3(d)} ~ 7(d) 
at d > 1.5. For D < D* (1 < d < 1.333...) this denomi- 
nator can be approximated by { 7 (d)}, because /3(d) = 
for d < 1.333. ...Table [I] gives a set of critical exponents, 
including (f3/v) e ff and (l/v) e ff, to be used below for 
scaling Monte Carlo data of ML 13 ^ vs leji 1 /^ for strips 
of various D and L = 500, 1000 calculated using periodic 
boundary conditions. 

Figures inland [3 give scaling plots of Monte Carlo data 
for strips of lenght L = 500, 1000 various widths D. 
Figure H3 shows results for D = 3 < D* ~ 6, taking 
|e| = |e~ 4 / T |, (one-dimensional scaling) and Figure Eire- 
suits for D = 20 > D* ~ 6, taking |e| = |T - T c | (two- 
dimensional scaling). 



III. CONCLUDING REMARKS 

In summary, our Monte Carlo data properly analyzed 
show that strips of dimensions LxD, with D <C L, 
are characterized by a susceptibility critical exponents 
"f e ff(D) blowing up at D = D* ~ 6, which appears to 
correspond to a dimensionality d* , 1.333... < d* < 1.5, 
which is well determined by the extrapolation of (3(d) and 
8 (d) given by the same linear expressions which de- 
scribe well the dimensionality dependence for d = 2, 3, 4. 
This lends empirical support to the theoretical contention 
that the renormalization transformations, which are not 
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TABLE I: Ising critical exponents for dimensionalities 



1< d < 2. 


deff 


Peff 


5 eff 


leff 






1.0 








1/2 


2 





1.1 








(0.71) 


(1.54) 





1.2 








(1.25) 


(0.96) 





1.3 








(5.00) 


(0.26) 





1.333... 








oo 








1.5 


0.031 





oo 








1.6 


0.050 


0.013 


3.700 


0.421 


0.021 


1.7 


0.069 


0.027 


2.516 


0.642 


0.044 


1.8 


0.087 


0.040 


2.100 


0.791 


0.069 


1.9 


0.106 


0.053 


1.887 


0.905 


0.096 



The numerical values are obtained by means of the linearly ex- 
trapolated expressions for -y(d) = (3(d)[S(d) — 1] , /3(d) = (|g - j), 
8-Hd) = (fjf - i), v(d) = y(d)/d and P(d)Md) = (j^). 
Values in parenthesis for D < D* are approximations obteined, 
extrapolating d e ff(D) from the D dependence at D > D* . 



only scale invariant but also conformally invariant trans- 
formations, imply fractional exponents for Ising lattices 
of geometry intermediate between one dimensional and 
two dimensional. 

It may be pointed out that exact information for the 
largest eigenvalues of the tranfer matrix of the system 
in small width strips (perhaps up to D < 20 or so with 
infinite length) can be obtained using f.i., the Lanczos 
algorithm, but, in principle, this does not result in useful 
information about the effective exponents, which are the 
ones which determine the observable behaviour of M(T) 
close to, but not arbitrarily close to the transition. 
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FIG. 2: Inverse susceptibility vs |e| = |e 4//T | for L x D strips 
with L — 1000 and D = 2, 3, 4, 5 showing evolution of ')eff(D) 
for D <D* ~ 6. 




FIG. 3: Inverse susceptibility vs |e| = \T — T c \ for L x D strips 
with L = 1000 and D = 12, 16, 20, 24; showing evolution of 
leff(D) for D>D* ~ 6. 
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FIG. 4: Effective dimensionality d of L x D strips as a function 
of D for _D > D* ~ 6 obteined by means of the relationship 
7 e //(D) = 7(d) = /3(d) [5(d) — 1] as given by means of linear 
extrapolations of /3(d) and 5~ 1 (d) from the known fractional 
values at d = 2 and d = 4 (see text). 
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FIG. 5: Plot of susceptibility critical exponent 7e//(-D) for 
L x D strips as a function of D for strips with various widths 
D < D* and D > D* ~ 6. The continuous curve for D > D* 
is given by J e ff(D) = G + yl(D - D*) -1 - S(D - D*) with 
G = 11/8, A = 1125/24 and B = 3/3200 in good agreemente 
with 7(d) = /3(d) [8(d) — 1] as determined by from /3(d) and 
8 (d) linearly extrapolated from the respective d = 2, d = 4 
values. 
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FIG. 6: Scaling plots of ML 13 '" vs \e\L 1/u for strips with 
D = 3 < D* taking |e| = \e~ 4/T \. 



